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SUMMARY 


A new approach to the design of multivariable control systems using 
the Multivariable Nyquist Array method has been developed. The technique 
utilizes a conjugate direction function minimization algorithm to achieve 
a diagonal dominant condition over the extended frequency range of the 
control system. The minimization is performed on the ratio of the moduli 
of the off-diagonal terms to the moduli of the diagonal terras of either the 
inverse or direct open loop transfer function matrix. 

In addition to the ability to achieve diagonal dominance with a mini- 
raiam of designer intervention/ several new feedback design concepts and 
evaluative measures are introduced. These include: 

1. Dominance control parameters for each control loop. 

2 . Compensator normalization to evaluate open loop 
conditions for alternative design configurations. 

3. An interaction index to determine . the degree and type 
of system interaction when all feedback loops are 
closed simultaneously. 

This new design capability has been implemented on an IBM 360/75 in a 
batch mode but- can be easily adapted to an interactive computer facility. 
The design method represents a significant contribution . to the design and 
analysis of multivariable control systems in the frequency domain and has 
been applied to the- Pratt and Whitney FI 00 'turbo fan engine with three 
inputs and three outputs. 
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SECTION 1. 


INTRODUCTION 

Over the past several decades, considerable effort has been expended 
in the development and synthesis of linear multivariable feedback control 
theory and its application to the design of multivariable control systems. 

Initially the. analysis of automatic control systems utilized time 
functions usually expressed in differential or integral form. Maxwell, 
more famous for his work in field theory, presented the first mathematical 
treatment of a control mechanism in 1868 CH* During the two decades 
preceding World War II, important contributions took place in aviation, 
electronics and circuit theory. For example, Nyquist's classic work in- 
1932 [2] on stability -of linear feedback systems was prompted not by sta- 
bility problems in control theory but by a desire to better understand the 
characteristics of certain communication networks. During World War II 
these concepts were rediscovered by control people and have since played an ' 
important role in the- control field . 

By 19M-5 the theory of linear servomechanisms and the fundamentals of 
mathematical modeling were well developed. The concept of steady state 
transfer functions had been introduced by Harris [3] and incorporated into 
the earlier work of Nyquist to further the understanding of the dynamic 
behavior and design of servomechanisms. This mathematical concept was popu- 
larized by Gardner and Barnes [4U with the introduction of transform 
calculus. By the end of the forties, the analysis and synthesis of linear 
continuous systems was basically limited to trial and error methods. 

Around 1950, Evans [5] introduced the Root Locus method, which for the 
first time provided a means for the direct sjmthesis of control systems. 
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These techniques have since been further developed and still represent one 
of the most useful synthesis techniques available for linear systems; 

t 

In the mid-1950’s, control engineering sees an unprecedented growth 
rate. Analog, digital and hybrid computers have reached high levels of 
perfection and more important, become universally available. The single 
input single output frequency domain techniques of Bode-, Nyquist and Evans- 
can now be coded for computer generated plots and the world of computer 
aided system design is created. 

With the computational speed and numerical accuracy now available and 
inspired by the work of Russian and American research teams, the field of 
control enters a new era referred to as "modern control theory" encompass- 
ing the general areas of optimal and adaptive control. Vector space methods 
provide the mathematical foundations for the time domain synthesis of multi 
input multi' output con-trol systems. In 1960, Kalman C6] provided a defini- 
tive trea-tment of the linear case with a quadratic cost function and showed 
that the optimal feedback control is determined by the unique positive semi 
definite solu-tion of the matrix Riccati differential equation. A special 
issue of the IEEE Transactions on Automatic Control Theory (1971) on the 
Linear Quadratic Gaussian (LQ6) problem reflects the tremendous breadth and 
depth of this field. 

The modern con-trol era is often characterized as the "algorithmic era". 
The system description, design goals and parameter constraints are in many 
cases manipulated and massaged until the problem format -fits a description 

t 

for which there is an algorithmic solution. The resulting controls are 
usually highly interactive, require full state feedback- and generally result 
in low integrity systems. Often the problem has been so narrowly defined 
that the final .control configuration is unique. 
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In general, fhe success of linear optimal control theory (and other 
pole placement algorithms), when viewed from the frequency domain, is that 
full state feedback ensures adequate gain and phase margins in each of 
the feedback loops. If all system states are not available for meas\ire- 
ment, severe penalties in terms of phase lag may be incurred- Techniques 
have therefore been developed to provide estimates of the inaccessible 
states [7,8]. This, in turn, leads to a dynamic feedback controller. 

In the late I960’s, it became clear to Rosenbrock, MacFarlane and 
others that the vector time response methods leading to the LQR problem 
and associated regulator solution methods were not the panacea long pro- 
mised. Optimal control design techniques although suitable for multi input 
multi output system analysis did not possess many of the design capabili- 
ties of the classical methods. By a suitable generalization of the 
frequency response methods, originally introduced by Bode, Nyquist and 
Weiner, to the multi input multi output system, new dimensions in the 
classical design concepts were created. 

An algebraic theory, based on Rosenbrock' s work, defines the struc- 
tural relationships in terras of which feedback systems may be manipulated 
into a variety of feedback forms. The generalization of Nyquist 's funda- 
mental criterion, the concept of integrity and Bode's sensitivity results 
were extended to vector forms. A survey of the major results in linear 
multivariable feedback theory from the vector frequency response viewpoint 
was outlined by MacFarlane [9-11]. Of those techniques presently available, 
the Inverse Nyquist Array [12] introduced by Rosenbrock and the Characteri- 
stic Locus Method [13], introduced by MacFarlane have surfaced as two of the 
most •• useful frequency domain design techniques for a wide range of practi- 
cal multivariable feedback systems. Both methods require a computer-aided 
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design facility with an interactive graphic display unit upon which the 
appropriate loci are computed and displayed. 

After an initial inquiry into the design capabilities, mathematical 
dependencies and computational requirements of the multivariable frequency 
domain techniques,' it becomes apparent that the basic principles and the 
governing philosophy of the Inverse Nyquist Array (INA) provide for the 
maximum utility of the single loop classical design theories. The INA 
extends Nyquist’ s stability criterion to inverse polar plots and multi in- 
put multi output systems. It provides the mechanism to reduce system 
interaction to a degree wherein each feedback loop can be independently 
designed. It utilizes the theorems of Gershgorln and Ostrowski to deli- 
neate the bounds of the eigenvalues of the transfer matrix at each fre- 
quency and thereby define the degree of both the open loop and closed loop 
interactive effects. These basic concepts and the underlying components 
associated with an INA design are further outlined in section 2. 

The single most detracting feature of the INA design philosophy is 
the unreasonahly-high degree of designer intervention to. secxare the condi- 
tion of ’’diagonal dominance". Since this condition is crucial to the 
success of an INA design, a principal objective of this research project 
was to develop an alternative method to search for the dominant condition. 
In section -3, of this report, a new algorithm utilizing a conjugate direc- 
tion' function minimization technique is presented. In fact, the algorithm 
is sufficiently versatile so as to be appropriate to the design of a. multi 
input multi output system using the Direct Nyquist Array (DNA) in addition 
to the INA. 

, In addition to the tremendous versatility and flexibility of the -pro- 
posed algorithm, several new design concepts are introduced in section 3. 
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These new concepts respond to the concerns pertaining to the level of 
-dominance required, to complete the design on a single loop basis, the de- 
gree 'and type of closed loop interaction which will result upon simultane- 
ous closure of all feedback loops when each loop has been designed inde- 
pendently and a method of comparing system compensators each of which 
produce the desired dominance condition. 

In section several 'examples are presented to demonstrate the compu- 
tational efficiency and effectiveness of the proposed technique. The 
principal example of this section is the analysis and design for the Pratt 
and -Whitney FIDO- turbofan engine at sea level static conditions. 
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SECTION 2 


MULTIVARIABLE NYQUIST ARRAY 

The Multivariable Nyquist Array (MNA) design method is herein pro- 
posed as the union of two mutually exclusive design techniques : the 

Inverse Nyquist Array (INA) and the Direct Nyquist Array (DNA) methods. 
’Both methods have identical design objectives and are founded upon a com- 
mon mathematical structure. The methods are mutually exclusive in the 
sense that the INA utilizes inverse polar plots while the DNA -uses direct 
polar plots. The principal point of departure is the use and interpreta- 
tion of the multivariable Nyquist stability criterion in achieving the 
final system design. 

The fundamental objective of the MNA design methods is to decrease 
system interaction to such an extent that the closed loop system design 
.problem reduces .to a set of independent single loop design problems. 
Although simply stated, the actual reduction procedure proposed by 
Rosenbrock [12,14] requires a high degree of designer intervention and is 
fundamentally a trial and error procedure. . The algorithm developed in 
subsequent sections of this report considerably reduces this designer 
dependency thus making the MNA design method a more viable design tool. 

Historically, the' first attempt to eliminate system Interaction was 
proposed by Boksenbom and Hood [15], Their procedure was to completely 
decouple system input output pairs through the appropriate design of pre 
and post compensator matrices. The resulting compensator forms are 
necessarily complicated and in many cases unstable and/or physically un- 
realizable. It did, however, provide for single loop closure on a 
completely non- interactive basis. Further attempts at system decoupling 
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are reported in Cl6] and [17], 

The Multivariable Nyquist Array method .adopts a considerably more 
sophisticated viewpoint in an attempt to achieve a similar design condi- 
tion. The MNA recognizes that the extreme degree of decoupling theory is 
not really necessary to establish the desired design conditions. Compen- 
sator designs should be stable and realizable and" preferably as simple as 
possible for ease of implementation. It further recognizes that some 
degree of system interaction may actually be desirable in the event of 
sensor or actuator failures. 

In the remaining parts of this section, the mathematical foundations 
of the' MNA are.- briefly introduced followed by an outline of the INA and 
DNA design methods. The section concludes with a discussion of the advan' 
tages and limitations of MNA desl^ philosophy as originally proposed by 
Rosenbrock. 

A. Mathematical Perspectives 

In Figure 2.1, 6(s) is an mxm transfer ma-trix representing the coup- 
ling of the m inputs to the m outputs. The pre and post compensator 



Figure' 2.1 Multivariable System Configuration. 


7 






matrices, K(s) and L(s), respectively, are each of dimension mxm. The 
feedback gain matrix F(s) is assumed to be diagonal of similar dimensions, 
•Clearly, if ■ 

Q(s) = L(s)G(s)K(s) 

is diagonal, loop closure may proceed on an individual loop basis with a 
guarantee of zero system interaction. It is this premise upon which the 
MNA is based. However, the adherence to s'trict diagonalization is relaxed 
with the substitution and exploitation of the concept of diagonal dominant 
matrices . 

Definition A matrix is said to be diagonal dominant when the moduli of 
' its diagonal elements are greater than the sum of the moduli 
' • of the . corresponding off-diagonal elements, taken by row or 
by column. 

That is, if Z is a mxm complex ma'trix then 

a. Z is diagonally row dominant if 

m 

jz. .j- Z \z..\ > o for all i = l,2,...,m 
j=l 

b. Z is diagonally colTimn 'dominant if 

m ■ 

lz..|- E |z.,| >o for all i = l,2..,,Tn. 
j'=l . 

It is important to note that to satisfy the definition Z must be entirely 
row dominant or entirely column dominant. The definition does not provide 
for a mixture of row and colvimn dominance for different diagonal elements. 

Quite apart from the concept of dominance is a theorem by Gershgorin 
which states that all eigenvalues of a complex matrix Z are located in the 
union of the circular discs defined by 

]a - z.. I 1 r i = l,2,...,m 
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with -center at and radii given by 


or 


m 

r • = Z I z. . 

1 j,=l 

iT^j 


m 

I \z..\ 

Di' 


(by rows) 


(by columns) 


i=l 


The latter radii follow from the fact that the eigenvalues of Z are equal 
to the eigenvalues of Z. 

Now let Z be a function of the complex variable s, i.e. 



z^iiCs) 


• • • 



22i<=) 




Z(s) = 

• 

• 

* 


* 

m 



2.1 


At each value of s on a specified contour D in the s-plane, Z(s) is a com- 
plex matrix and the preceding definition and theorems apply at each and 
every point on D. Thus the eigenvalues of Z(s) are functions of s' and the 
concept of diagonal dominance can be reformulated as 


m 


11 


lz..(s)|- Z jz.,(s)| > o (row dominance) 
j=l 


2.2 


and 


m 


z. .(s) - Z Z. .(s) > o 


11 


j=l 




(column dominance)- 


2.3 


Let D be a large contour in the complex s-plane consisting of the 
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imaginary axis from s = -jR to s = +jR together with a semicircle of 
ra<Jius R in the right half plane as indicated in Figure 2,2. As s tra- 
verses the D contour in a clockwise direction, z..(s) will generate a 



Figure 2.2 the D contour in the s-plane 

curve for i = l,2,...,m' in the complex plane. From the application of 
Gershgorin’ s theorem at each point on D, a band of circles centered about 
will similarly be generated as in Figure 2.3. If a separate figure is 
constructed for each loci, then the collection of figures represent a set 


of "fuzzy” Nyquist (or inverse Nyquist) diagrams, one for each input-output 
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Clearly, the eigenvalues of Z(s) are captured within the union of 
the Gershgorin bands since the eigenvalues at each point lie within the 
corresponding union of discs. If the matrix Z(s) is diagonal for all s 
on D, then the Gershgorin band reduces to zero width and the are the 
characteristic Ipci of Z(s). Thus the width of the Gershgorin band pro- 
vides a qualitative measure of the departure from the diagonal condition. 
In a control system setting this departure would reflect the degree of 
open loop system interaction. 

A further interpretation of Gershgorin ’s theorem states that the 

eigenvalues of.Z(s.) lip in the intersection of the union of Gershgorin ■ 

bands by. row and by column. This follws immediately from the eigenvalues 
T 

of Z (s) being equal to those of Z(s). Thus if the intersection excludes 
the origin then the determinant of Z(s) cannot be zero, i.e., if Z(s) is 
diagonal dominant then det | Z(s) ] ?^o. With this corollary the following 
theorem results Cl4]. 

Theorem 2.1 Let Z(s) he an mxm rational matrix and let D be a closed^ ele- 
mentary contou r hay ing, on it no pole of z^^(s), i = 1,2,. ..,m. Let there 
exist e>:b such that for each s on D either 

m 

[z..(s),|“ E |z..'(s)j > e i = 1,2,..'. ,m 
. j=l . 

or 

m 

[z. .(s) j- S |z. -(s) I > e i = 1,2,-. ..,m 

j=i 

Let z^^(s) map D into i = l,2,...,ra, and let det|z(s)| map C into 

r^. Let encircle the origin times, and let encircle the origin 

N times . Then 
z 
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m 


N 


i=l 


N. 

1 


This theorem provides for a number of stability conditions for the control 
problem depending upon which matrices are dominant. Further considerations 
are given in the individual discussions of the MNA methods. 

In review, - it is the Gershgorin theorem bounding the eigenvalues of a 
complex matrix and the concept of ’diagonal dominance which provide the 
necessary foundations for the HNA methods. Theorem 2.1 will form the basis 
for closed loop stability considerations. Final system design will then 
proceed on an individual loop basis. 


B. Inverse Nyquist Array 

To in’troduce the Inverse Nyquist Array, consider Figure 2.1 to repre- 
sent a block diagram for a single input-single output feedback control 
system. For this system, the open loop transfer function is 

Q(s) = L(s)G(s)K(s) ■ 2,4 

Let F(s) be a constant scalar feedback gain with closed loop system trans- 
fer function H( s ) , given by 


H(s) 


Q(s) 

l+Q(s)F 


2.5 


The INA method uses inverse relationships for a variety of reasons, one 
of which is the simplification of equation 2.5, i.e,. 


H 


^ l+Q(s)F 

ts; - -^33- 


= Q“ (s) + F 


2 . 6 " 


Here the inverse .closed loop transfer function is simply, a linear transla- 


tion, in the complex • plane of the inverse open loop transfer function. 

For notational convenience and to avoid confusion at later stages of the 
development the following definitions are made 
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H(s) 

Q(s) 

K(s) 

' Thus, (2.6) becomes 

-H(s) 


Q"^(s) 

K“^(s) 

A 

Q(s) + F 


2.7 


2.8 


Now let Q(s) map the D contour* (in the s-plane) into which encircles 

• A A • A 

the origin N times clockwise. Also let H(s) in (2.8) map D into F^^ 

A 

which encircles the origin times clockwise.. Here D is the usual 
Nyquist contour in the s-plane and is sufficiently large to include all 


finite poles and zeros of Q(s) and H(s) in the closed right half plane. 

For the system of equation (2.8), N, is the number of times F encir- 

M 

cles' the point 


(~ F, o) 

and the following statement of the Nyquist Stability Criterion for the INA 
results : 

Theorem 2 . 2 Let the open loop system, Q(s), have poles in. the closed 
right half plane. Then the closed loop system is as 3 miptotically stable if 
and only if 


N - N, = p 
q h o 

A ' A 

Figure 2. '4 represents a typical INA polot of F^. Each crossing of F^^ 
by the critical point represents the entry of a pole of the closed loop 
system into the right half plane if the crossing is from "left to right". 
A "right to left" crossing represents the removal of a pole from the 
right half plane. All directional crossings are relative to an observer 

A 

on F in the direction of increasing frequency. 
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Figure 2.4 Typical INA plot of F^. 

Some important observations regarding the INA plot of Figure 2.4 can now 
be made: 

1. The INA begins at zero frequency and terminates at infinite 
frequency. This characteristic provides for an immediate display of 
low frequency information essential to most design problems. 

2. With F=:0, -the origin (0,0) represents the open loop situation. 

Thus if the system is open loop stable, then any gain in the vicinity 
of the origin which does not cross represents a stable closed loop 
system operation. 

3. The closed loop inverse transfer function H(s) is given by the same 
diagram with the origin shifted to C-F,0). 

4. Only positive frequency is plotted with the negative frequency 

range inferred. Thus if the feedback gain were such that the critical ■ 

point was at point A, then two poles are about to cross into the right 

half plane at the gain crossover frequency o) . . 

A 
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5. The ratio OA to OB is the gain margin of the- closed loop system 
with gain F. 

6. The phase- margin of the closed loop system -is y when a gain of F 
is introduced into the feedback path. 

' A . /S 

7. The steady state offset is given- by the ratio Q(0)/H(0) and is 
obtained from the INA plot directly as the ratio of OC to BC. As F 
increases, the offset decreases. 

8. The effective bandwidth of the closed loop system is given by the 
value of CO for which 

. lH(jw)|- = 1.414|F(jC0)| 

9. Fee(pDack compensation procedures for the INA are' conceptaually 
similar to the direct Nyquist array. For example, to provide more phase 
advance a phase lead compensator could be introduced into the feedback 
path. In Figure 2.4, this compensator would shift the INA plot up and 
away from the origin with gain crossover occ\u?ring at point E. 

Clearly, the Inverse Nyquist ^ray for single -input single-output systems 
is at least as versatile a design tool as the direct Nyquist method. Many 
authors contend that the simplification of the closed loop transfer func- 
tion coupled with the low frequency profile make the INA diagram a more 
useful design mechanism for single loop design.' A more detailed study of 
inverse polar plots is contained in Rosenbrock [14] and Raven Cl8]. 

The design ‘of feedback control units for multi input multi- output systems 
using' the Inverse Nyquist Array is relatively straightforward once the con- 
dition of diagonal dominance has- been achieved. Fundamentally, this condi- 
tion suggests that system interaction has been reduced to such an extent 
that each control loop can be closed separately and independently from the 
remaining loops using .single loop theory . 
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Let Figure 2.1 represent a feedback control system with m inputs and 
m outputs. The system transfer matrix G(s) is mxm with each element 
g..(s) consisting of a numerator and denominator polynomial representing 

J 

the transmittance between input j and output i. In general, G(s) or 

/s 

,G(s) will not be dominant over the frequency range of interest and there- 
fore must be modified to confo™ to the requirements of the multivariable 
array methods. 

For the INA method, the pre and post compensator matrices (each mxm) 
must be selected such that 

Q(s) = K(s)G(s)L(s) 2.9 

is diagonally dominant over the D contour. Using the definitions 'of 
(2.2) and (2.3), Q(s) is diagonal dominant if 

a. E |q. . (s) j / |q. . (s) I < 1 (row dominant) 2.10 

» _ in n 

1=1 

i^j ■ • 

or- 

b. . E |q. . (s) I / I q. . (s ) I < 1 (column dominant) 2.11 

for all i = l-,2,...,m. From Gershgorin’s Theorem, all eigenvalues of 
Q(s) are captured within the union of the Gershgorin bands centered about 
with radius as the sum of, the moduli of the off diagonal terms taken 
by row or by column. 

Assume, for the moment, that Q(s) is diagonal - dominant by rows over the 
D contour as in Figure 2.5 for m=2. If the system is open loop stable then 
the closed loop system will be guaranteed asymptotically stable for all 
feedback gains on the real axis in the vicinity of the origin and bounded by 
the Gershgorin bands. This stability criterion (Theorem 2.1) can now be 
restated as follows; 
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i 


a . • Row 1 


b. Row 2 


Figure 2.5 Q(s) diagonal dominant with m=2. 


Theorem 2.3 Let each of the Gershgorin bands based on the diagonal ele- 

/\ A. 

ments q^^(s) of Q(s) exclude the origin and the point (-f^,0). Let these 


bands encircle the origin times and encircle the point (-f^,0), 
times. Then the closed -loop system is asymptotically stable if and. only if 


m 

2 

i=l 



m 

Z 

i=l 


N. 


hi 


= P, 


2.12 


where is the number of open loop poles of Q(s) in the right half plane. 

Here the Gershgorin bands are defined by the radii of (2.10) or (2.11). 

Theorem 2.3 is stated in its most useful form for application purposes 

since K. . and N, . are evaluated from the same set of Gershgorin bands, 
qi hx 

The theorem could be stated in a more general form wherein new bands would 
be recalculated for each set of f^, 

i • 

Note that theorem 2.3 specifically states band encirclements and not just 
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encirclements of the T . . Thus the theorems and corollaries for the INA 

3 . 

apply only for feedback gain values located outside the Gershgorin bands. 

The stability theorems simply do not apply for any gain values located 
within a band and ho inferences regarding stability or instability can be 
made. 

For the system in Figure 2.5, each feedback loop could be independently 
closed using the inner most envelope of the Gershgorin band for gain margin 
and phase margin assessment. In general, the Gershgorin band provides a 
conservative estimate of the stable gain space and is most useful as a 
preliminary design tool. 

To provide further insight into the INA design mechanism, assume' Q(s) 
diagonal dominant and all feedback gains f^ chosen in accordance with 
Theorem 2.3. Let hi,^(s) be the transmittance from input i to output i when 
all feedback loops are closed i.e. f^J^o, i=l,2,...,m. Define r^(s)- as the 
transmittance from input i to output i when all feedback loops, except the 
ith loop, are closed, i.e. f^=0. From standard feedback relationships 

hii(s) = r^Cs)/(l+r^(s)fjj^) . 2.13 

hi'^(s) = r^(s) + f^ 2.14 

Hence, to complete a .set of single loop designs by opening one feedback 
path at a time, it is the quantity r^(s) which governs the system behavior 
and not Rosenbrock exploited this relationship to demonstrate that 

r^(s) is located within the Gershgorin band for all stabilizing feedback 
gains. He further demonstrated that when all gains except f^ are specified, 

A 

the transmittances ^^(s) for i=l,2,..., m are located within a narrower set 
of bands. This new set of bands is based upon a theorem by Ostrowski and 
are appropriately labeled the Ostrowski bands. For the INA method, the 
Ostrowski bands are always located within the Gershgorin bands. 
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Theorem 2.4 (Osirrowski [19]) Let D be the closed Nyquist contour in the 

s-plane and let Z(s) be .row dominant on D with no pole of z^^(s) 

i=lj2,,..,m on B. Then if is a point on D, Z(s^) has an inverse 

Z(s ) such that for i=l,2, . , ..,m 

}z..(s ) - z.. ^(s )| < <i. d, (s ) < d.(s ) 

. ‘ 11 O 11 o' ^1 1 o 1 o 


when 


m 


d.(s^) = E |z^j(s^)| i = l,2,...-,m 

ji^i 

are the Gershgorih radii and are a set of "shrinking factors" defined 
by 

d, (S q) 

4.(s„.) = max I- ■ V — vr 

.10 . z. .(s ) 

i ' ii o ' 

j?*i 

A yv. A 

For the INA,Z(s) = H(s) = Q(s) + F, hence 

|h^^ ^(s) - < <|>j^(s)d^(s) 

for row dominance, and 

Ih^/^Cs) - Cf^ + q^^(s)-)|< <{>^(s)d^(s) 

for column dominance. The shrinking factors thus become 


d.(s) 

IfTflTTSTT 


and 


ij).(s) = max I - — f — ^ 

1 I • iCi « • \ S J I 

' D 99 ‘ 

Using (2.14), ,(2.15) and (2.16) become 

|r.^ - q^^(s)| < (j)^(s)d^(s) < d^(s) 

for row dominant 'H(s) and 

/ 

Ir^^ ^ - q^^(s)| < <fi^(s')d^(s) < d^(s) 


2.15 


2,16 


2‘.17 


2.18 


2.19 


2.20 
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for column dominant H(s). 

The OstrowsJci bands serve two useful functions. First they locate 
the transmittance s r^ within a narrower set of bands when all loops except 
the ith-are closed. As the f^ j?^i are varied, r^ will also vary so that 
the ith Ostrowski band depends upon f^ jp^i. Second, they provide a more 
accurate measure of phase and gain margin for the ith control loop and 
thus reduce the problem to the design of f^. An extremely important feature 
of theorem 2.4 is that once the f^ l^j are specified to obtain the Ostrowski 
band for loop j , the band for loop j will continue to shrink when the 
f^ ij^j increase.’ That is, if the feedback control in loops i?^j are in- 
creased to improve the control' in the respective loops, the feedback design 
for loop j is unaffected. , . - • 

The Inverse Nyquist Array method consist of the following fundamental 
operations : 

1. Design K(s) and L(s) such that Q(s) is dominant by row or column 

using ( 2 . 10 3 of ( 2 . 11 ) . • ' 

2. Plot the Gershgorin Bands for each control loop. 

3. Evaluate stability for the diagonal dominant' Q(s) using theorem 


2.3. 

4. Finalize the design using "the Ostrowski Bands and single loop 
con'trol theory. 

Clearly, the application of the INA method is predicated upon the 
ability of the system designer to 'achieve an adequate degree of system 
dominance. The degree of dominance attainable is primarily governed by two 
factors : 


1. The structural sophistication and realization ofpre and post com- 
pensator forms;. 
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2. The ability of the designer to manipulate the compensator para- 
meters to achieve, maintain or improve dominance in each row (or 

-s 

column of Q(s). 

Current methods available for compensator design require a high degree 
of designer interaction and are therefore best suited for interactive 
computer facilities. Often success in achieving dominance is based upon 
the experience and intuition of the system designer in the application of 
these trial and error methods. In addition the methods are restricted, 
for the most part, to the design of constant compensators although attempts 
at more sophisticated structural forms has been reported [143. 

One of the first methods suggested beyond a total trial and error 
approach was to diagonalize G(s) at s=0 [12] by setting. 

K = LG(0)' 

If G(0) is non-singular, Q(s) will be dominant near the origin of the com- 
plex plane but not necessarily over the extended dynamic frequency range 
of the system. ■ In either event, however, no general guidelines exist for 
improving dominance in any given row (or column) of Q(s) since any modi- 

A. A 

fication of the coefficients of K (or L) translate to- vector addition when 
forming 

A. At A* A, 

Q(s) = KG(s)L 

Similar conditions exist for high frequency dlagonalization using 

K = LG(oo) 

Therefore dlagonalization at either end of the frequency spectrum provides, 
at best, a starting point for the trial and error approach. 

The pseudodiagonalization method developed by Hawkins [20] and genera- 
lized by Rosenbrock [14] formulates the dominance objective as an eigen- 
value-eigenvector problem. The method is best suited for constant compen- 
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sator design but could be used for polynomial forms of K(s). Briefly, 

A 

pseudodiagonalization determines the set of row elements of K (for row 

A 

dominance) or column elements of L (for column dominance) which "most 
nearly" diagonalizes Q(s) at specified frequencies. For the eigenvalue- 
eigenvector problem this translates to minimizing the weighted sum of 

A 

squares of the off diagonal elements of Q(s) at the N specified frequenc- 
ies, 

N m ^ • 2 

Min 2 T { S lq--(jw )1‘ } i = l,2,...,m 2.21 

r=l j=l 

where m is the dimension of G(s) and are designer specified weighting 
factors. -The minimization in (2.21) is subject to either a constraint on 
the elements (for row dominance) 


m 

1 



1 


2.22 


or a constraint on the diagonal element 

lq^.j(j(fl)| = 1 (0 specified 2.23 

If the constraint in (2.22) is used, q^^ may vanish since it does not 

appear in the problem formulation. If N=1 in (2.21) and constraint (2.23) 

/ . 
is used, problems analogous to the case when (u=o may result. The additional 

feature here, however, is that the eigenvalue- eigenvector problem may be 

resolved at any frequency within the range of interest. The solution at the 

selected frequency is then tested at other frequencies for dominance. For 

N>^2, the eigenvalue-eigenvector problem resulting from the minimization - 

of (2.21) subject to (2.22) (or (2.23)) further complicates the dominance 

issue since oj in (2.23), N values of oj and N values for y in (2,21) must 

r r 

be selected a priori. If either pseudodiagonalization scheme yields domi- 
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nance over the frequency spectrum then the designer may resort to trial 
and error methods to improve the degree of dominance. However, if domi- 
nance does not prevail, conclusions concerning the non-existence of con- 
stant compensators may be erroneous . This would be the case, for example, 
when constant compensators exist which yield dominance over the D contour 
in the sense of (2.10) or (2.11) but did not satisfy (2.21) for any combi- 
nation of to , X and N. 
r’ r . 


C. Direct Nyquist Array 

The direct Nyquist Array (DNA) method for single-input single-output 
systems is well established in the control literature. Fimdamentally, it 
is a polar plot of the open loop transfer function as s traverses the 
familiar Nyquist contour. 

If Figure 2. '2 represents a single input single output feedback system, 
then the open loop transfer function is 

Q(s) = L(s)G(s)K(s) 2.24 

with the closed loop transfer function as 


H(s) = 


Q(s) 


2,25 


l+Q(s)F(s) 

Let (|)(s) be the characteristic polynomial of H(s), i.e,, 

(Ks) = 1 + Q(s)F(s) 2.26 

and let <Ji(s) map D into a closed contour F. as s traverses D in a clock- 
wise direction. From the Nyquist stability criterion, if encircles the 
origin times clockwise then the closed loop system is asymptotically 
stable if and only if 


N, = 0 2.27 

<i> 

An equivalent but more convenient form of the Nyquist criterion is 


available using (2.24). Let Q(s) map D into the closed contottr F. and let 

9 
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F(s) be a constant different from zero. Taking the usual precautions to 
insure that D encloses every finite pole and zero of (j)(s) in the right 
half plane, (2.26) can be rearranged to 


ii>(s) = F + Q(s) 2.28 

and the following theorem results Cl4] 

Theorem 2 . 5 Let F be constant and- let the open loop system have poles 
in the closed right half plane. Let Q(s) map D into f making N clock- 
wise encirclements of the point (- -^,0). Then the closed loop system of 
Figure 2‘.1 is as 5 nnptotioally stable if and only if 

N = -p 2.29 

. q ° 

In theorem 2.5, (2.28) could be written as 
N = Z - p 

q 0*^0 

where Z^ is the number of finite zeros of \p(s) in the closed right half 
plane. However, the system is stable if and only if Z^=0 thus (2.29) re- 
sults. A typical Nyquist polar plot is indicated in Figure 2.6 for, posi- 
tive frequencies. The. closed loop system design can now be completed using 


gain margin, phase margin, etc. as the criteria and - jr as the critical 



Figure 2.6 Nyquist polar plot. 
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point'. Further design considerations for the single input single output 
case can be found in [18] and [21]. 

The extension of the- above design concept and in. particular theorem 
2.5 to multi input multi output systems will require the diagonal dominant 
condition imposed on Q(s).' Therefore, let Q(s) be diagonal dominant, i.e.. 

Row dominant 2.30 

3. “■ 1 , 2 , . . . ,TR 


Column dominant , 2.31 

X — l,2,..«,m 

for all s on D. Let F ^ (or equivalently F) be a diagonal constant matrix, 
■fchen the off diagonal elements of 'f'(s), 

f(s) - f“^+Q(s)- 2.32- 


m 

kii(s)| > iq^.(s) 


j=l 


or 


m 


j=l 


are equal to the off diagonal elements of Q(s). The diagonal elements 
of 2.32 are 

■ ' ' • ' 2.33 


T|>ii(s) = f^^ + q^^.Cs) i = 1,2,. 
--1 


,m 


Now let f.. + q. .(s) map D into a closed curve and let q..(s) 

11 ^11 ^ f+q. . ^11 

^11 

-1 

map D into F^^^. From Gershgorin's Theorem, all eigenvalues of F + Q(s) 

are captured within the union of the Gershgorin hands centered about 

F^ ^ with radius equal to the sum of the moduli of the off diagonal 
“^ii 

terms taken by row or by column. But the Gershgorin radii- for F ^ + Q(s) 
are equal to the radii for Q(s). Therefore exclusion of the origin by a 
band centered about F^ ^ is eq^ui valent to the exclusion of the point 

C-f.^,0) by a band centered about F . The following theorem summarizes 

1 q-. 


this thought [141. 
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Theorem 2.6 Let each of the Gershgorin bands centered about exclude 

the point Let these bands encircle the point 

“1 

times, i=l,2,...,m. Then the closed loop system is 
asymptotically stable if and only if 


m 

S N. = -p 
. , 1 o 

1=1 


2.34 


Theorem 2.6 provides the necessary criteria to complete the DNA de- 
sign on a single loop basis. To improve the design capability through the 
reduction of the Gershgorin band, a theorem analogous to the Ostrowski 
theorem for the INA is available Cl4]. 

Theorem 2.7 -Let F be a diagonal constant matrix, and let F ^+Q(s) be 
dominant on D. Let h^^(s) represent the transfer function from input i 
to output i when all feedback loops except the ith loop are closed. Then 
for each s on D, 


Iq^^(s) - h^(s)l < if'^(s)d^(s) < d^(s) 

for row dominance" and 

kii(s) - h^(s)( < (}>j^(s)d^(s) < d^(s) 

for colTimn dominance where 

<ji^' = max 


2.35 


2.36 


2.37 


.j. jf..+q..(s)| 

1 

0.(s) = max — 


2.38 


and the d-(s) and d. (s) are the appropriate Gershgorin radii. 

J J 

Extreme care must be exercised in using theorem 2.7 in a DNA applica- 

i 

tion. This is evident, for example, in the row dominant conditions of 
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(2.35) and (2.37). Here it is observed that the shrinking factors in 
(2.37) decrease as the decrease, as they must since 

. k^^(s)| = |h^(s)[ 2.39 

in the open loop condition. Therefore, in contrast to the INA, the 
Ostrowski bands increase as the feedback gains in the remaining .loops in 
crease. This condition may cause some difficulty if the feedback design 
for loop j is based upon the set {f ij^ j } and some of the f^ increase as 
a result of later design modifications. The increase in value of some of 
the f^ will cause the Ostrowski band for loop j to increase which in some 
cases may be sufficiently large so as to encircle the design point in loop 
j.. If this situation prevails, the entire DNA design is voided, since the 
critical design point in every loop must remain outside the corresponding 
Ostrowski band for Theorem 2,6 to be valid. Hence the feedback design in 
each loop must be. reevaluated whenever significant positive increments are 
made in the remaining feedback loops. 

The above situation is cited by Rosenbrock [14] as the single most 
determining factor for choosing the INA method over the DNA method. How- 
ever, the graphical interpretation of theorem 2.7 is not quite as grim as 
Rosenbrock may suggest. The following reasons are cited: 

1. For any specified set of feedback gains [f^, i?^j], the transfer 
function from input j to 'Output j is contained within the jth 
Gershgorin band (Theorem 2.6). Thus if the Gershgorin bands are 
sufficiently narrow then single loop closure can be completed without 
invoking theorem 2.7. 

■ 2, Theorem 2,7 can be used in precisely the same sense as theorem 2,4 
was used in an INA design, if the feedback gains {f£,ii^j} are un- 
, usually large (and outside the respective Gershgorin Bands), then as 
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the are reduced under single loop closure, the Ostrowski bands 
in theorem 2,7 will shrink. 

3. Examination of (2.37) (or (2.38)) suggests that if all ij^j 

1 ' 

are infinite, then the largest value of <|)(s) (or (ji (s)) is given 
by 


max 

j 


d^(s) 


2.40 


But this expression is simply the maximum degree of dominance 
attained in the remaining Gershgorin bands. Therefore each Gersh- 
gorin band may be immediately reduced by the appropriate factor in 
,(2.40). This important concept is used in the algorithm proposed in 
section 3 of this report. 

Apart from the graphical interpretation and use of Ostrovski's Theorem, 
the DNA and the INA methods are similar in concept. For every theorem, 
definition and design concept pertaining to the INA method there exists an 
analogous , theorem, definition and design concept for the DNA method. The ' 
methods' are mutually exclusive in the sense that a design initiated using 
the INA method cannot be completed using DNA design concepts and vice versa. 
The main point of departure lies in the graphical interpretation of the 
fundamental theorems, wherein the DNA method utilizes polar plots to 
interpret the design objectives and the INA utilizes the inverse polar plot. 

Clearly, the DNA method (as well as the INA method) is critically 
dependent upon the ability of the system designer to achieve a diagonal 
dominant condition. All methods appropriate to the INA method to achieve 
a dominant condition are also suitable to the DNA design' method. 

The DNA method can be summarized as follows; 
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1. Design K(s) and L(s) for diagonal dominance of Q(s) by row or by 
column. 

2. Plot the Gershgorin bands. 

3. Evaluate system stability using theorem 2.6. 

4. Con^lete the design using single loop control theory within the 
limitations of theorem 2 . 1 . 


C. Discussion 

The multivariable Nyquist array is a very useful and versatile design 
technique for multi input multi output systems. The principal feature of 
the MNA is the utilization of the mathematical foundations of complex 
matrices to reduce system interaction to the extent that each feedback 
loop can be independently designed. The Gershgorin and Ostrowski theorems 
are easy to apply and the final closed loop design concepts are well docu- 
mented in the control literature. 

i 

The utility of the MA design philosophy is totally dependent upon the 
ability of the system designer to achieve an adequate degree of system 
dominance. It is therefore in the interests of the designer to have avail- 
able new methods of achieving the dominant condition and thus provide 
greater flexibility in the use of the MNA design methods. 

Ideally, any new dominance method sho^lld be sufficiently flexible to 
address the following issues: 

1. Design of .constant compensators for either row or column domi- 
nance . 

2. Design of frequency dependent compensators with fixed structural 
. form for row or column dominance. 

3. Eliminate the need for designer oriented trial and error’ methods 

to improve dominance. 
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4. Provide a measure for the comparison of two or more compensator 
I pairs each yielding a dominant condition. 

5. Utilize the degree of dominance in the remaining rows or columns 

to improve the design in a particular row or column. 

6. Be appropriate for both the INA and the ,DNA methods. 

This list highlights several major concerns regarding the use of the 
MNA which has not. been treated in the literature to date. .The concerns 
are primarily directed toward the degree of dominance required in a 
particular design, the utilization of dominance to further improve the 
design base for any given row or column element and the' degree and type of 
system interaction resulting from the simultaneous closure of all feedback 
loops. 

A new technique to generate diagonal dominant compensator pairs is 
described in the. following sections. The method addresses each of the 
concerns cited above and is suitable for implementation in either an 
Interactive or batch computer mode. 
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SECTION 3 


NEW DOMINANCE ABGORITHM 

In this section, a new dominance algorithm for the Multivariable 
Nyquist Array is developed. The algorithm utilizes the explicit defini- 
tion of' diagonal dominance of section 2 to minimize a non-anal 3 rtic 
function of the pre and post compensator ma-trix parameters. A conjugate 
direction function minimization technique applied over the dynamic fre- 
quency range of interest is used to achieve the desired dominant condition. 

In its most general form, the new dominance algorithm is characterized 
by four specific phases: 

a. Parameter initialization 

b . Dominance evaluation 

c . Parameter optimization 

d. Design considerations 

The most unique features of the algorithm are the non-interactive nature of 
phases b and c and the increased degree of flexibility and designer inter- 
action in phase d. In addition, the algorithm is sufficiently general so 
as to be appropriate to the batch computer mode design of both the Inverse 
Nyquist Array and the Direct Nyquist Array. 

The algorithm is developed in the following manner: First, the con- 

cept of dominance is reviewed from a computational and structural viewpoint. 
The main computational unit, which performs phase b and phase c above, is 
then introduced followed by the program con-trol unit. Finally, closed loop 
system design concepts new to the MNA design philosophy are presented and 
discussed. 
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A. Dominance Observations ' 


The general concept of diagonal ‘dominance was presented in section 2. 
It is‘ now of interest to carefully examine the detailed structure of the 
dominance condition and to identify the points of similarity for row and 
colimih dominance determination. 

Consider a general complex matrix Z(s) to be m x m and the square 
matrices A and B to be constant. Let P(s) be the matrix product 

P(s) = A Z(s) B 3.1 

with diagonal dominance of P(s) defined as 
m 

^ ^ i— Ij 2 ^»..,m 3.2 

j=l ^ 

m 

for row dominance and 
m 

Z ip..(s)i/lp^^(s)l < 0. < 1 i=l,2,...,m 3.3 

j=l ^ 
ij^j 

for column dominance. In (3.2) and (3.3), 9^ represent a specified set of 
constants contained in the semi-open interval (0,l3 and will be defined in 
part D of this section. 

‘Using (3.2) as the definition of row dominance and assuming B to be 
specified, the following observations can be made: 

OBS 1 

Diagonal dominance for row i of P(s) is determined exclusively by the 
elements of the ith row of A^ 

OBS 2 

Dominance of row i of P(s) is unaffected by a scaling of row i of A. 
OBS 1 states that the elements of row i of the A matrix (a.. j=l,2, 
...,m) do not enter into the consideration of dominance for any other row 
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Hence, only m parameters in the A matrix need be determined for each 
row dominant condition. OBS 2 suggests that the row elements of A can be 
multiplied by a common factor with a guarantee of dominance preservation. 
This observation is used in part D to provide a mode of comparison of 
compensation forms and is fundamentally a scaling or normalization opera- 
tion. 

OBS 3 

Using (3.3) with B specified, column dominance of (3.1) may be ob- 
tained only if all elements of matrix A are manipulated simultaneously. 

This observation follows immediately upon noting that all elements of 
A, influences the behavior of every- column of P(s). As a result of 

OBS 3, the method of pseudodiagonalization cannot be used to search for 
column dominance when B is specified. 

The transpose of P(s) in (3.1) provides a similar set of conditions 
for the B matrix when A is specified. Specifically, let 

P^(s) = p\s) = B*^ z"^(s) A*^ = A^ Z^(s) B^ 3.4 

In (3.4), Aj^ = B*^, 2j^(s) = z"^(s) and = A*^. Hence, observations 1-3 
are now appropriate to (3'.4) and are summarized as follows 
OBS 4 

If A is specified in (3.1), then column dominance of P(s) may be 
attained via the manipulation of the elements of B wherein the elements of 
the jth column of B only Influence the dominance condition in the jth 
column of P(s). 

OBS 5 

If A is specified in (3.1), then row dominance of P(s) may be obtained 

2 ^ 
only if the m elements of B can be simultaneously manipulated. 
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The above observations are clarified upon examination of P(s) in 
(3.1) when m=2 

^11^22^12^®^ ^12^22^22^®^ 

^ 21 ^ 12 ^ 11 ^®^ '*’ ^ 22 ^ 12 ^ 21 ^®^ 

+ a2j^i>22^12^®^ +- a22b22Z2l^®^ 

3.5 


-If the elements of both A and B are unspecified, no obvious pattern 

for Obtaining diagonal dominance of P(s) is identifiable beyond the simul- 

2 

taneous manipulation of the 2m parameters. In addition, the elements of 
B are, in a sense, competing with the elements of A in the attempt to 
secure dominance of P(s). 

In view of OBS 6, it is apparent that designer intervention may be 
required to establish a hierarchial structure to the dominance evaluation 
procedure. This particular situation has not been automated in the pro- 
posed algorithm and remains a subject for future research. Several guide- 
lines are presented in the discussion of the examples of section 4. 

In the conditions imposed above, it was implicitly assumed that con- 
stant compensators were the desired form of system compensation. This 
would certainly be true from an implementation perspective and is there- 
fore the initial assumption in any MNA design attempt. However, in some 
applications, constant compensator matrices may not yield the dominant 
condition or satisfactorily meet the closed loop design specifications. 


P(s) = 


^ 11 ^ 11 ^ 11 ^®^ ■*' ^ 12 ^ 11 ^ 21 ^®^ 


+ ^12^21^22^®^ 


® 21 ^ 11 ^ 11 ^®^ ^ ^ 22 ^ 11 ^ 21 ^®^ 


® 21 ^ 21 ^ 12 ^®^ ^ 22 ^ 21 ^ 21 ^®^ 


Further examination of (3.5) yields 
OBS 6 
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The alternative is therefore dynamic compensation. 

The observations made above for constant compensators can be extended 
to include the design of dynamic compensators. To illustrate, consider 
the objective to be row dominance ofP(s) in (3,1) with a specified post- 
compensator matrix B(s). A(s) in (3.1) can be generalized to 


A(s) = 


®22^22^®"' ^2m^2m^®' 






a f (s) 
Iran mm 


3.6 


where fj^^(s) are specified by the system designer and the a^^ are to be 
determined for dominance of P(s). With respect to the a^^’s, OBS 1-6 are 
retained. A special case of (3.6) is the pre- compensator form 


A(s) = [A^ + Aj^/s] 


3.7 


OBS 7 


In an IKA design using dynamic compensation, it may be desirable from 
an implementation viewpoint to structurally define the compensator in the 
inverse domain. 

This observation is explored more fully in part D, but serves as a. 
reminder that if A(s-) in (3.7) is the precompensator form for the IMA, then 
A ^(s) could have poles and zeros in the open right half plane. This 
situation arises from the general observation that 


a.. f..(s) f [a.. f..(s)] 
1] 1]'“ " 13 13 


-1 


3.8 


and thus could be avoided with the precompensator form 

A(s) = [A^ + A,/s:"^ 

o 1 


3.9 
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Here OBS 3 would be imposed for both row and column dominance attempts. 

V?ith this insight into the structureLl and computational form of 
dominance evaluation, it is quite evident that (3.1) could represent 
either the INA or the DNA open loop design format. Thus a computational 
unit devoted to the search for diagonal-dorainant-forming. compensator 
parameters can be developed. 

B. Generalized Optimization Unit 

The main unit of the dominance algorithm is constructed in a gener- 
alized setting and is founded upon the observations in part A. 

Let C represent a vector formed by consolidating the unknown compen- 
sator parameters into a single array. In some applications, such as row 
dominance of P(s) with B .specified', C will be an m vector representing 

the unspecified values in the jth row of A. For other . situations , re- 

2 

fleeted in OBS 6, C may be of dimension 2m . Ass^lme further that for any 
given C vector, P(s) can be properly reconstructed and evaluated for all 
s. The main unit then performs the following functions: 

1. Accepts proper coding to identi:^ the MNA design, form and select 

the desired performance measure for dominance evaluation. 

2. Adjust the elements of C to minimize the performance measure 

selected using a conjugate direction function minimization technique. 

f 

The performance measure selected is dependent upon the form of domi- 
nance desired. Fundamentally, there are two specific forms from which 
the selection is made: 

m 

a. J. (C,d. ) = Max Z |p * -(to) [ / |p . . (u) j 3.10 

j=l 

i?fj 
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3,11 


m 

b. J(C, 6 • ) = Max Max C S |p . . (w) [ / |P . . (oj) j - 6 • ] 

■ ^ (oeSJ i j=l ^ 

The measure in (3.10) is selected whenever the unspecified coefficients 
can be subdivided into mutually exclusive sets, ^en this situation 
occurs m calls to the main unit will transpire, one call for each value 
of i. Alternatively, the measure in (3.11) is used whenever the dominance 
seeking coefficients must be simultaneously manipulated. Details are 
postponed -to the following subsections. 

The success of any numerical optimization technique in locating the 

1 

extrema of a function of many variables is highly dependent upon the • shape 
of the contours of the function to be extremized and the convergence pro- 
perties of the optimization technique employed. Clearly, these concerns 
are imbedded in the successful evaluation of the compensator parameters 
in (3.10) or (3'ill). 

Examination of (3.11) reflects the interesting computational and 
numerical aspects of the performance measure. Each of the p..(o))*s are 
complex functions of the compensator parameters from which the definition 
of dominant matrices is composed. At each frequency, the rows (or 
columns) are scanned to identify the largest ratio. This array is then 
scanned (over w) to determine its maximum value. If the largest ratio is 
less 'than 9^^, then the trial elements generating P(s) yield dominance for all 
rows (or columns) over the frequency range of concern. If the ratio is 
greater than 0^, then the elements of C must be adjusted to create the desired 
condition. 

The majority of numerical optimization techniques reg^uire some form 
of localized gradient or second variational calculation. If the perform- 
ance-measure is a well defined anal 3 rtical function, then a gradient 
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dependent algorithm is often the most efficient route to pursue. This is 
the premise upon which the method of pseudodiagonalization is based. 
However, gradient calculations for (3.10) or (3.11) could result in 
significant numerical difficiiities due, in part, to the extremization over 
(D. For this reason a numerical optimization technique which does not 
explicitly depend on localized gradient calculations is preferred. 

The method selected for implementation is the Zangwill-Powell [22,233 
optimization technique. This method is known to be effective when sharp 
ridges and narrow valleys are present in the performance contours and is 
suitable for application to problems with a large set of variables to be 
optimized [22,24,253. Other methods which do not require gradient calcu- 
lations and therefore might be appropriate for the optimization of (3.10) 
or (3.11) are those of Swann [263 (an extension of Rosenbrock [273), 

Smith [28 3 and Wood [293* Of the optimization methods which do not re- 
quire derivatives, Fletcher’s study [303 suggests that Powell’s method may 
be computationally the most efficient. 

Fxmdamentally, the optimization unit of the dominance algorithm per- 
forms as follows: 

' 1. Upon receipt of an initial guess for the C vector, P(s) is 
evaluated to determine the ratio of the Gershgorin radii to the 

i 

moduli of the corresponding diagonal element in (3.10) or (3.11) 
for all Msfl. 

2. The dominance ratio is appropriately scanned over w to identify 
the largest ratio. The numerical values of J(C,9^) at subsequent 
evaluations are used by the conjugate direction minimization 
algorithm to adjust the components of C. 

3. Step 2 is repeated until either the desired degree of dominance 
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is achieved or internal checks within the optimization method indi- 
cate no further improvement in J(C’,0^) is likely with successive 
adjustments of C. 

This optimization unit forms the nucleus of the dominant seeking al- 
gorithm. Access to and control of this Tinit is the function of the program 
supervisor.. 


C, Program Supervisor 

The program supervisor performs the task of accepting information from 
the system designer and properly coding the optimization unit and the de- 
sign unit (see next subsection) to perform the requested design. The 
supervisor will accept the following data for each design attempt: 

1. INA or DNA design? 

2.. System dimension 

3. Frequency range 

4. Frequency increment (equal spacing) 

5 . • Row or column dominance? 

6. Identify, fixed compensator 

i. Precorapensator 

ii. Postcompensator 

iii. None 

iv. Both 

7. Numerator and denominator coefficients for each element of G(s) 

8. Dominance control parameters, 6^., i = l,2,..,,m 

9 . Plot options 

i. No plot 

ii. Row or column elements superimposed by rows only 
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iii; Gershgorin Land only 
XV. (ii) and Ciii)- only 
V, Gershgorin and Ostrowski bands 
vi» Ostrowski only 

10, Diagonal feedback gain elements for Ostrowski bands 

11 . Precompensator specifications 

i. Dynamic or constant? 

ii. If dynamic, specify coefficients of f..(s) 

iii . Initialize coefficients 

a. Set to identity matrix 

b. Diagonalize at 

c. Elements to be read in 
.. 12. Post compensator specifications 

i. Dynamic or constant? 

ii. If d 3 mamic, specify coefficients pf f^^(s) 

iii. Initialize coefficients 

a. Set to identity matrix 

b. Diagonalize at 

c. Elements to be read in 
13, Advance design control parameters 

With the above information the MNA design is completely specified and 

• -9 

a search for dominance using the optimization unit may be implemented. In 
the following subsection, design concepts new to the MNA design philosophy 
are introduced. -Section 4’ illustrates the use of these concepts in appli- 
cations . 
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D. New Design Concepts 


The Inverse Nyqiiist Array as conceived by Rosenbrock and implemented 
on a PDP-10 digital computer by Munro at the University of Manchester was 
intended to be used in an interactive computer mode. In this mode, the 
system designer is an integral part of the computational and evaluative 
phase of dominance determination. As a result of this high degree of 
designer intervention, final design considerations are based upon trial and 
error methods which, in most cases, are not systematized. It is the intent 
of this subsection to Introduce new computer aided design techniques to the 
multivariable Nyquist array. 

1. Dominance Control Parameters 

Theoretically, any degree of system dominance is sufficient for the 
application of the'MNA design method. The degree of system dominance may 
range from the marginally dominant condition (performance meas\ire less than 
but near unity.) to the decoupled condition (performance measure near zero). 
Correspondingly, if the level of system interaction is interpreted in terms 
of the width of “the Gershgorin bands then there is a direct correlation be- 
tween system interaction and dominance. 

From a practical viewpoint, the degree and type of system interaction 
is an important design consideration in the selection of" input-output pairs 
and corresponding compensator structures. For this reason it may be 
desirable to reduce open loop system interaction before the feedback loops 
are closed. This is accomplished through a specification of the "domi- 
nance control parameters", 6^ i=l,...,m in the system performance measures 
of (3,2) and (3.3) (or alternatively (3.10) and (3.11)). 

In application, the system designer will specify each 6^ for i = 1,2, 


. . . ,m where 




3.12 


0 < 6 ^ £ 1 

The unspecified parameters in the compensator matrices are then- adjusted 
by the optimization unit in an attempt to meet this degree of dominance. 

The 9^ (selected in accordance with (3.12)) is fundamentally a request to 
the optimization unit to make the largest Gershgorin radius smaller than 
100 9^ percent of the corresponding diagonal element. If the optimization 
unit -can satis:^ the 0^ request in each row (or column), then the prescribed 
degree of dominance has been achieved. 

An interesting observation regarding Ostrowski’s theorem can be made 
when. the dominance control specifications are satisfieds Recall that the 
shrinking factors in Ostrowski’ s theorem are 

V 

' d (s.) 

(INA) (j).(s) = max 7-^ 3-13 - 

X ...■ f.+z.-(s) • 


and 


(DNA) <J).(s) = max 

j?*i 


dj(s) 


3 


Letfj = 0 in an INA design and f^ = oo in a DNA design, then 


<I>£(S') 


d.(s) 

... I z . . ( s ) I 
3^1 ‘ 33 ' 


max 


3.15 


Define 


9. 

1 


Max 

s 


<l>£(s) 


Max 

s 


Max 

3i^i 


d^.(s) 


hence 


3.16 


^^(s) £ 9^ for all i = l,2,...,m. 

But the are nothing other than the largest degree of dominance obtained 
for j?^i, therefore 
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a. 

1 


3.17 


= max 0 . 

■ 1 

3?i 

Thus for any MNA design, the Gershgorin band may be immediately reduced by 
the corresponding 9^ factor. 

2. Interaction Index C31] 

Once a suitable degree of dominance has been obtained, each feedback 
control loop may be independently closed using single loop classical con- 
trol theory as outlined in section 2 of this report. Application of 
Ostrowski’s theorem can then be used to further reduce the Gershgorin 
bands for each design loop. The corresponding Ostrowski bands can thus be 
used as a conservative estimate of the design parameters to improve the 
overall closed loop system design. 

It is important to re-emphasize that each loop is designed indepen- 
dently. Furthermore, the width of the final set of Ostrowski bands, in a 
broad sense, reflect the degree of closed loop system interaction in the 
finalized design. However, no information regarding the type of interaction 
is available from the Ostrowski plots. For this information, the "Inter- 
action Index" developed by Davison [31] is employed. 

Briefly, the interaction index assumes that m linear time invariant 
proportional feedback control loops have been independently designed and 
are separately applied to the system. Davison considers the question of 
how much interaction will occur when all m control loops are to be applied 
simultaneously to the system. 

. Consider the linear time invariant system 

x(t) = A x(t) + B v(t) 3.18 

y(t) = C x(t) 3.19 

In an MNA format, the control law y(t) is 
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• y(t) = - K F L C x(t) + yO(t) 3.20 

where K and L are the dominance producing compensators, the feedback 
gains have been determined from the Ostrowski plots and y'’(t) are the 
input disturbances to the system. 

Under the assumption that m control laws have been found so that the 
resultant systems are stable and the jth controller satisfactorily con- 
trols output y^ , the following index of performance is chosen 


J. = max j y. (t) dt j = l,2,...,m 

^ X o ^ 
o 


3.21 


X X = 1 

o o 


Now,' define the interaction index as 


I. = 
J 


J. - J. 

_J 1 

J. 

3 . 


3 — 1 , 2 ,. ..,m 


3.22 


Here is the value of (3.21) resulting from, .the application of y. only 

and Jj is the value of (3.21) when all loops have been closed simultan- 
eously. The index in (3.22) thus provides a measure of the relative change 
in the control of when all feedback loops are simultaneously closed 

compared with the jth control law applied independently. 

The problem defined by (3.18), (3.19), (3.20) and (3.21) can be re- 
formulated in a Lyapunov function: 


J . - max X Q.x = A (Q.) 

] o ^3 o max 3 

.0 


3.23 


X X = 1 

o . o 


where A (Q.) is ,the largest eigenvalue of Q. and Q. is the solution to 
max 3 ' 33 


CA - BKF.LC] Q. + Q. [A - BKF.LC] = " C , C.' 

3 3 3 3 3 3 


3.24 
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The Cj in (3.24) are the corresponding rows of matrix C. The matrix 

in (3,24) has zeros in all positions, except for fjj» which is the feed- 

& 

back gain for loop j . To evaluate is diagonal with each diagonal 

element set to the corresponding design gain. 

Using (3,23), the interaction index becomes 


I. = 


(Q. ) 

max 

X — (qTT 

max 3 


- 1 


3.25 


and will lie between the bounds 


- 1 < I . < «> 
1 


3.26 


a. - 1 < I. < 0 
3 

A 

This situation occurs when the performance index for is less than 
and is construed to be a favorable form of system interaction. That is, 
the control of output is in^foved when, all feedback loops are simul- 
taneously closed when compared to the control of y^ using only input . 

The ultimate limit of -1 is obtained when a high degree of system inter- 

A 

action exists such that J. is near zero’, i.e. 

1 


J. < <’J. 
3 3 


3.27 


b. I. = 0 

3 

This situation occurs when 
A 

J. = J. 3.28 

3 3 

and thus implies that the closure of the remaining feedback loops has no 
effect upon the control of y^. This case corresponds • to a decoupled condi- 
tion. 

C. I. > 0 
3 

Here loop closure has a deleterious effect upon the control of y. when 
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compared to the single loop closure of Vj . Clearly, this form of system 

interaction is undesirable when I. becomes large. As I. tends toward 
» 3 3 

infinity , the interaction becomes more severe and implies a tendency • 
toward instability (1^ = “). 

The interaction index thus becomes an important evaluative measure 
for closed loop system design and can be applied directly to the MNA. 


3. Comparison of Compensator Designs 
In the • application of the Multivariable Nyquist Array design method, 
alternate structural forms and different initial guesses for the unspeci- 
fied compensator parameters may yield (after optimization) different 
compensator designs, each of which satisfy the dominant conditions. It 
is therefore of interest to have a means for the comparison of compensators 
which may have been generated by different design forms. 

In particular, consider an INA design with L specified and K con- 
strained to be constant. For row dominance of Q(s), each row of K is 
mutually exclusive of the elements of the remaining rows. Thus the ith 
row of K may be normalized about the diagonal parameter. This procedure 

A, 

applied to alternate designs for K will provide a direct comparison of 
dominant compensator forms and therefore may be used to determine the 
suitability and/or superiority of the alternate pre- compensator design 
forms. 

As an illustration, consider an mth order G(s) with L = I. Assume 
that m initial guesses for K yield m different normalized compensators 

A 

r = l,2,...,m each yielding dominance for all rows. of the corresponding 

A 

Qj,(s). For each row the best design is selected and the corresponding row 

A A ' 

of identified. The composite K matrix is then used to finalize the 
design. 
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This procediire could also be used for other MNA design forms and thus 
provide the desired mode for comparison. 

For dissimilar structural forms each yielding an acceptable design, 
the Interaction Index discussed above could be used as the mode for 
comparison. Here the comparison of row versus column designs or designs 
using different input-output pairs could be made. Thus any final design 
decision could ultimately be based upon the level of system interaction and 
actuator and/or sensor failure accomodation. 

E. Discussion 

In this section, a new algorithm to obtain the diagonal dominant con- 
dition for the multivariable Nyquist array has been developed. The algo- 
rithm is compatible to both the INA and DNA design philosophy with either 
constant or dynamic compensation. 

Fundamentally, the algorithm is based upon the observations and 
generalization of the concept of diagonal dominance. The characterization 
of dominance as a function minimization problem provides for a tremendous 
degree of design flexibility and eliminates a substantial portion of the 
trial. and error aspects of previously used methods. In addition, the pro- 
posed algorithm may find dominance conditions when other techniques are 
inappropriate or have previously failed. This is certainly true whenever 
a design attempt is made for which pseudodiagonalization cannot be used, 
i.e., INA design for column dominance with specified L matrix. 

Several advanced design concepts new to the MNA design philosophy have 
been introduced. These concepts are based upon the availability of a fast 
and efficient method to generate the dominance condition. The dominance 
control parameters are used in an attempt to secure a specified level of 
dominance in each control loop. They are used in the generalized optimi- 



zation unit and correspond to the minimization of the largest Gershgorin 
radii within the ’ frequency range specified. Since the frequency range is 
determined by the system designer it could vary from a single isolated 
point to a fixed interval in the frequency spectrum to the' entire spec- 
trum. The total impact of this degree of flexibility has yet to be 
realized. 

Utilization of the interaction index and compensator normalization 
methods suggest a means by which the system designer can evaluate a pro- 
posed closed loop design in addition to providing a comparison of competing 
designs. These techniques could also be used in the assessment of dynamic 

feedback components and thus provide a quantitative measure of closed loop 

t 

system response characteristics and interaction. 

In the -next section, the results obtained by the dominance algorithm 
are compared with previously reported applications. In addition, an analysis 
is performed on the FlOO turbofan engine using an INA format. In each 
case the primary goal is to either verify the reported results and suggest 
new alternatives or to simply obtain the dominance condition with little 
effort devoted to attaining a final design. 
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SECTION 4 


•APPLICATIONS 

The principal objective of this section is to demonstrate the effi- 
ciency of the algorithm in section 3 and to stibstantiate the suitability 
of the MNA design philosophy of section 2 as a viable closed loop design 
alternative for air breathing propulsion systems. In each of the appli- 
cations considered' herein, the principal concern has been to achieve the 
diagonal dominant condition with little attention devoted' to any final 
design considerations. 

The section is subdivided into three subsections: 

A. Previous applications of MNA 

B. F-lOO Turbofan Engine 

C. Discussion 

In part A, the algorithm is- applied to several test cases with new and 

interesting results to be reported. Part B is the application of the MNA 

to a- linearized operating point model of the Pratt and Whitney FlOO series 

2 engine at sea level static conditions. This represents the first attempt 

at an MNA design for a sixteenth order FlOO model. Part C provides an analysis 

of the application areas and suggests new dimensions for closed loop design of 
air breathing' propulsion systems. 

A. Previous applications of MNA 

At the present time, the control literature contains approximately seven 
or eight reported applications of the Inverse Nyquist Array. The design 
applications using the INA have been primarily conducted at the University 
of Manchester Institute of Science and Technology (UMISTO under the auspices 
of Prof. Rosenbrock and performed by his collegues (Munro, Rutherford, etc".) 
and students. Little activity regarding the INA has been reported outside 
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UMIST' due in part to the dependency of the IM on computer-aided graphic 
facilities. Since the algorithm of section three reduces this dependency, 
all previously reported design applications served as test cases for the 
new design algorithm. -In every case, the algorithm operating in the 
appropriate mode for evaluation purposes either confirmed the reported 
design compensator ' forms and their indicated numerical values, suggested 
alternate dominant forms or modified the reported compensator values to 
extend the dominant condition to a larger frequency range. 

The results of the trial cases were generally obtained in one pass 
through the algorithm using K = I or K = LG(0) as the initial guess with 
the specified frequency range stibdivided into N equally spaced points. 
Depending upon the dynamic frequency range of interest, N could be 
selected as any integer value between one and one thousand. On the 
average, each test case required 100 CPU seconds on a batch mode IBM 360/ 
75 to achieve a dominant condition and plot the indicated Gershgorin and/ 
or Ostrowski bands, 

This is in direct contrast to the typical two weeks to nine man- 
months of effort required to achieve dominance using the UMIST computer 
aided design suite. 

In each of the following test cases, the dominance algorithm of sec- 
tion 3 provided, new and interesting results regarding the specific appli- 
cations. 

1. Rosenbrock [14] - Boiler Furnace 

In this application, the objective is the control of a boiler 
furnace with four inputs and four outputs as represented by 
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• 1 

.7 

'.3 

. 2 

1+4 s 

i+5s 

1+5S . 

1+5S 

.6 

_ 1 

.4 

.35 

1+5S 

1+4S 

1+5 s 

1+5S 


.4 

1 

.6 

1+5S 

1+5 s 

1+4S 

1+5 s 

.2 

.3 

.7 

.,1 

l+5s 

l+5s 

l+5s 

l+4s 


For this G( s ) , ’ Rosenbrock has determined that the system is diagonal domi- 
nant in the INA sense (G(s))- when 


K = L = I 4.,2 

It is further demonstrated that the dominance condition is improved if. 

>v 

L = I arid 


K = G(0) = 


1 

,6 

,35 

.2 


.7 

1 

.4 

.3 


.3 

.4 

1 

.7 


.2 

.35 

,6 

1 


4.3 


Using pseudo-diagonalization at ta = 0.9 for all rows with L = I, the 
following precompensator is obtained 


1.469' 

79,44 

-.148 

.050 

-.654 

1.814 

-.249 

-.229 

-.229. 

. -.249 

1.818 

-.654 

.050 

-.148 

-.944 

1.469 


The effect of C4.3) is to reduce the Gershgorin band to the point where the 
system is essentially non-interacting. 
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Using the algorithm of section 3, the following runs were made in 
one submission of the computer program; 

a. K=Ij L=I; INA; No optimization 

b. Same as (a) except K=G(0) 

c. Same as (a) except K was set equal to Rosenbrock's solution 

d. K=I; L=I; INA; optimize with 6^=0, for i=1.2.3,4; set all 
feedback gains to zero and plot Gershgorin and Ostrowski ' 
Bands 

e. Same as (d) except K=G(0 ) 

f. Same as (d) except K = Rosenbrock Solution 

In runs d and e, the dominance control parameters were set to zero to 
determine, the set the precompensator values which will yield a condition 
of maximum dominance . . Runs ■ a through c were made to verify the reported 
conditions in [143 - 

'for each run, system dominance was obtained. The following table 
provides a comparison of the degree of dominance achieved in each case: 

Table 4.1 


Degree of 
Dominance 
Achieved 

• Run 
a 

Run 

b 

Run 

c 

Run 

d 

■Run 

e 

Run- 

f 


.88041 

.24899 

.25402 

.13118 

;13134 

.13117 


. 81731 

.28357 

.24924 

; 13841 

.14184 

.13865 


. 81731 

.28357 

.24924 

.13814 

.13846 

.13835 


.88041 

.24899 

.25402 

.13684 

.13340 

.13471 


A comparison of columns a through c verify the reported conditions namely 
that , pseudodiagonalization (Run c) reduced the degree of system interaction 
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to a significantly lower level- However, an overview of Table 4.1 inci- 

cates that the proposed algorithm was able to further improve the degree 
/ 

of dominance by approximately fifty percent (runs d and e). The final 
set of precompensator values are indicated below: 


Run d: INA, K=I- L=I,. 0.=O.O desired 




.91320 

.58078 

.25031 

.16693 




.49444 

.91399 

.33017 

.28856 








4.-5 



.28259 

.32363 

.90150 

.48944 




.16965 

.24860 

.58459 

.90534 


Run e: 

INA„ K=G(0), L=I, 

9^=0. 0 desired 

- 



1.10040- 

. 70007 

. 30000 

.20000 


K = 
e 


■ . 54478, 

1.01650 

.37337 

.32247 

4.6 



.34419 

.39708 

1.10420 

.‘59983 




. 19154 

.28148 

.65139 

1.0143 


Run f: 

•INA, K = 

Rosenbrock solution, L= 

=1, 9^ = 0.0 desired 

% 


.,84900 

.53952 

.23280 

.15500. 


y\ 

*^f 

- 

-.43153 

.79200 

.28384 

.24560' 

4.7 


. 

.23436 

.26755 

.74215 

.40167 




16510 

. 24714 

.59242 

.92996 


It is 

clear from 

table 4.1 

that any one 

of the precompensators corre- 

spending to 

runs d, e 

or f would be adequate 

for this application. However 

in other applications 

it may be 

necessary to 

form a composite compensator 


based upon the results obtained from different initial guesses for the 
unknown parameters. This is caused primarily by the inability of the opti- 
mization technique to detect and correct for the presence of local minima 
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during the optimization phase of the program. 


For this case. Table 4.1 would suggest the utilization of row 1 from 

in (4.7), rows 2 and 3 from K, in (4.5) and row 4 from K in (‘4.6) to 
f ’ d e - 

yield the composite form 


.84900 

.53952. 

.23280 

.15500 

.49444 

..91399- 

.33017 

.28856 

.28259 

.32363 

.90150 

.48944 

.19154 

.28148 

.65139 

1.01430 


with corresponding dominance levels as 

.13117 
' .13841 

'-e- = 

. .13814 

-.13684 


4.9 


From Ostrowski's theorem, the transmittance from input i to output i 
when all feedback loops are closed (except for loop i), can be located 

I 

within a narrower set of bands within the Gershgorin bands using the 
shrinking, factors 


<}>j^(s) 


jf.+q. ..(s)j 


d.(s-) 


4.10 


where d^Cs) are the Gershgorin radii, 
qjj(s) the diagonal element in Q(s). 
are bounded by 


fj the feedback gain in loop j and 


For all the shrinking factors 


0 < Max 
3?^i 


dj(s) 


3?^1 


6. a. 

1 , = 1 


where 0^ are the corresponding dominance levels. In this application, the 

3 . are 
1 
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.13841 


g _ .13814 

.13841 
.13841 

which implies . that the Gershgorin band may be immediately reduced by 
approximately 86.2 percent. The Ostrowski band in each loop will shrink 
further when the' feedback gains are increased. 

The res.ults presented above were obtained in one batch mode submis- 
sion of the computer program with a total expenditure of 30 CPU seconds 
including program compile time. 


2. Munro [32] - Aircraft Autostabilization 
In this application, the control unit for a two input two output 
model of a delta-winged aircraft is developed using the Inverse Nyquist 
Array method with constant pre and post compensation matrices. Here the 
authors restricted the post- compensator to the form 


II 

0 

11 

or 

II 

1 

0 . 


° • 1 



0 

^22 


Using the Lj^ form, row dominance was obtained using pseudo diagonalization 
with the following results for 0 < o) < 4.0 


K 


- 12.08 -34.19 

- '4.39 2.31 


4.12 


and 



1 0 

0 6.67 


4.13 
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It should be noted that pseudodiagonalization does not assist in the 
selection of the elements of the postcompensator matrix when row domi- 
nance of Q(s) is bo be obtained via the manipulation of the K elements. 
Thus to obtain (4.12) and (4.13) requires a hierarchial approach to para- 
meter selection. The results indicated above correspond to a set of domi- 
nance control parameters of C»86, .393. 

With the algorithm of section 3, the following six cases were pro- 
grammed for one batch mode submission: 


a. 


K = I 



1, 5, 10 


b. K = LG(0) 



1, 5, 10 


In each case, the Lj^ form was selected for comparison with (4.13). Iden- 
tical results would have been achieved if the L form had been used with 

3i 




The following results were obtained with 0^ = 0.0: 


Table 4.2 


Initial 

A 

4 . 

K 

1 

5 

10 

I 

■ Dominance 

Dominance 

Dominance 


Row 2 only 

Row 2 only 

Row 1 only 

LG(0) 

Dominance 

Dominance 

Dominance 


Row ■ 2 only 

Row 2 only 

Row 2 only 


From Table 4.2, it is clear that the desired condition is obtained from 
using a composite pre compensator formed using the appropriate 

A, 

dominance producing row of K as follows 
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1.0 


-.8388(10^) 


a. K. ... , = I 
initxal 


fxnal 


6- = [1.4244 .57973] 

3 


b. -K. , = LG(0) K . , 

xnxtxal fxnal 


= [.55472 1.8911] 


-1.8616 .99938 

-1.4625 -3.0661 
-.37956 -.62166 


c. Composite precompensator 




fxnal 



-1.4625' -3.0661 

-1.8616 .99938 

[.55472 .57973] 


L 


1 

b 


0 

10 


This application serves to. demonstrate that for some systems, efforts 
to achieve the dominance condition for Q(s) may be dependent upon the 
initial starting guess for K. "To explore this further, a parametric study 
was performed' on 9 ^^ for the starting conditions of Table' 4.2. The results 
of this study are reflected in figures 4.1 and 4.2. For each value of ^^ 2 * 
the dominance control parameters were set to zero and the precompensator 
values optimized. The level of dominance achieved for each case reflects 
a local minimum in the performance index and is indicated in the figures 
for each row of Q(s). 

Prom figures 4,1 and 4.2 it is clear that a composite precompensator 
matrix yielding the- desired compensator values can be obtained for any 
^22 "the range 

5-9. £ ^22 
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No attempt was made to identify the upper limit for £^2 since, this 
value can be deduced from the Gershgorin plots. 

Using in (4.13) the results of Munro are confirmed using the 
procedure indicated above and may be improved to the dominance levels 
of C .8317 .3867 1 with 

-14.59 -30.59 

- 4'. 303 2.31 

A A 

A review of the printout corresponding to K = L = I with no optimiza- 
tion reflects, the greatest departure from dominance occuring in both rows 
near to = 1» Selecting = 10 initial guess for the precom- 

pensator as 



^initial 


R fG{jl)] 
0 


dominance was obtained (via' the algorithm) for both rows simultaneously 


with 


K 


final 


-.13171 -.27622 

-1.0447 .56083 


and dominance levels of 

0 = t. 55.475 .57973] 

This result demonstrates that alternate methods of selecting initial 
parameter values based upon the systems characteristics are available. 


3. Munro [33] - Automotive Gas Turbine 
In this application of the'-INA a two input two output transfer 
matrix is used to describe the dynamic characteristics of an automotive • 
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.8065 -i- 2.64 -(15S+1.42) 

d^+1.15S+.202 s^+12.’8S^+13.6s+2.36 

I 

G(s) = 

■ 3..95s^it2;l2srfr4.9 7.14s^+25.8s+9.35 

s-^+9 . 155^+9 . 395+1 . 62 5^20 . 8s\li6 . 45^ 

+111.65+188 

gas tTorbine engine over the frequency range 0 _< to £ 25. From [33] ' 
column diagonal dominance is obtained with L = I and 

. 1.15 -.5175 

K = 

1.30 .415 

A 

Column dominance for Q(s) with a fixed post compensator requires 
continuous monitoring of both column dominance indices when any parameter 

A 

in K is adjusted. This form of the algorithm is outlined" in section 3 and 

A 

the above results were confirmed with an initial guess for K as the 
identity matrix for L = I and 0 ^ oj ^ 25. 

To examine the effectiveness of the- algorithm in section 3, the fol- 
lowing input data was provided to the algorithm: 

a. Row dominance desired 

b. L = I specified 

c. K = I as initial guess 

d. 0 ^ to ^ 25 subdivided into two hundred equally spaced points 

The following row dominant results were obtained within 30 CPU seconds on 
the IBM 360/75: - 

1 0 
0 1 
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2.9730 

-4.2414 

K = - 




12.507 

1.0247 

with 




tl 

CD 

.38007 


®2- ^ 

.42557 


representing the' maximum degree of system dominance in the respective rows. 
With this degree of dominance, the Gershgorin bands can be immediately 
reduced by 57. 443% and 61;993%, respectively, with zero feedback gains. 

Recalling the definition of dominance from section 2, it is clear that 
if the system is both row and column dominant at a particiilar frequency, 
then the smaller radius can be used for each row or column element in 
determining the appropriate Gershgorin radius. Examination of the computer' 
listing corresponding to the row dominant conditions above indicates that 
qil(s) is both row and column dominant over the range 0 _< oj < .5 and 
4.0 < 0 ) ^ 25, and q 22 (s) is both row and column dominant over 0 £ m ^ 25. 
This condition implies that the smaller radius can be used for both elements 
everywhere except • 5 m 4 . 0 . The remarks above concerning band reduction 
will still apply to the new radii. 

The Gershgorin and Ostrpwski bands for the row dominant configuration 
possess the same shape and form as those reported in. [33.] for column domi- 
nance . 


4. Sain [34]- - Turbofan Engine 

In [34], the two input two output transfer matrix representing the 
d 3 mamic characteristics of a turbofan engine is presented as 
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Although the Gershgorin bands are indicated for a column dominant 
Q(s) in [34], the corresponding compensator matrices were not provided. 
However, the authors suggest that column dominance was obtained during 

I 

a search for row dominance using pseudodiagonalization about u = 1. Using 
this method to achieve dominance they conclude that ’’though our examples 
are . . . introductory, they do serve to' show that typical jet engine models 
do not yield trivial dominance questions". 

With the algorithm of section 3 and G(s) above, row dominance was 
obtained for 

L\ = 

1 ■ 0 

K = 

0 1 ■ 

The corresponding dominance levels are 



• .0 = C. 848051 .406645] 

which can be reduced to . 

[.75329 .36597] 


using 
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0 


1 


L = 


1.2755 -.016011 

-.86064 ,77641 

Setting the feedback gains to zero, the Ostrowski factors are 
3 = [.36597 .75329] 

wherein the appropriate Gershgorin band can be immediately reduced. A plot 
of the corresponding Ostrowski band indicates an essentially decoupled 
system as all' lines are coincident over most of the frequency range con- 
sidered. 

For this application, six hundred points were used over 0 _< o) 200. 
Beyond 200 radians little dynamic activity occurs and thus was not used in ' 
the dominance evaluations. -However the above compensator matrices retain 
the dominant condition for w ^ 200. 

5 , Other Cases 

The dominance algorithm of section 3 has been used to examine the form 
and level of dominance in numerous examples in Rosenbrock's text and the 
current literature. In each case, the algorithm confirmed the reported 
results and improved the level of dominance using the dominance control 
parameters. In cases where the system was not dominant beyond a specified 
frequency, the algorithm adjusted the compensator values, to secure domi- 
nance over the entire spectrum. 

B. F-lOO Turbo fan Engine 

The engine tinder consideration is a Pratt and Whitney FlOO-PW-lOO after- 
burning turbofan. The FIDO is a low' bypass ratio , twin spool 
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axial flow engine with the following components; 

1. ■ Three stage fan driven by a two stage turbine 

2. Ten stage compressor driven by an air cooled two 
stage turbine 

3. Main burner with an annular chamber 

4. Variable area exhaust nozzle 

Using the non-linear dynamic simulation of the FIDO engine and the 
offset derivative method, a set of linear dynamic equations in state 
variable form for each of the thirty seven operating points is reported 
in C35il. For this study the sea level static (SLS) intermediate point was 
selected. This is in correspondence with zero Mach number, zero altitude 

V 

and a power level angle of 83°. 

The linear model at SLS intermediate is a sixteenth order system with 
the following state variables: 

= Fan Speed 
= Compressor Speed 
= Compressor Discharge Pressure 
= Interturbine Volume Pressure 
= Augmentor Pressiu?e 
= Fan Inside Diameter Temperature 
= Duct Temperature 
= Coti^ressor Discharge Temperature 
= Burlier Exit Fast Response Temperature - 
'= Burner- Exit Slow Response Temperature 
= Burner Exit Total Temperature 
= Fan Turbine Inlet Fast Response Temperature 
= Fan Turbine =lEnlet Slow Response Temperature 


^5 

^6 

^7 

^8 

^9 

^10 

^11 

^12 

^13 
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= Fan Turbine Exit Temperature 
Xj^g = Duct Exit Temperature, T^g^ 

Xj^g = Duct Exit Temperature, 

The engine inputs and outputs used for this study are: 

a. Two Input - Two Output Model 

Inputs: Uj^ = Main Bin?ner Fuel Flow 

U 2 = Nozzle Jet Area 

Outputs: y^^ = Fan Speed 

y '2 = Compressor Speed 

b. Three Input - Three Output Model 

Inputs: = Main Bilrner Fuel Flow 

U 2 ' = Nozzle Jet Area 

Ug = Inlet Guide Vane Position 

Outputs: y^ = Fan Speed 

y 2 = Compressor Speed 

yg = Augmentor Pressure 

The A and B matrices corresponding to the above models using 
X = -A X + B y 
y = C X 

are contained on Page 65 and 66 of C35], Application of Danielevski * s 
method for computing G(s) yields the following set of transfer functions: 

IS 14 1 

G^j^(s) = C- .0if57s t 54.789s + 45.(10 )s 

+ .1384(l0®)s^^ + .2211(0^°)s^^ 

T*? 10 1 U d 

+ ..20558(10 )s + .11944(10 )s 

+ .4534(10^^)3® + .11516(10^'^)s^ 
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+ .197(10^®) s® + .2246(10^^)3^ 

On u . . On 0 

+ .1657(10 )s + .75166(10 )s 

+ .1923(10^^)3^ + .2323(10^^)8 

4- .91419(10 )]/A(s) 

Gi2(s) =. [- 451. - ,3095(10^)3^'^ + .2062(10®)^^ 
1886(10^^)8^^ + .30163(10^^)8^^ 

+ . 2589 ( 10 ^^) 8 ^^ + . 1441 ( 10 ^^) 8 ^ 

-S- ,5592(10^^)3® + .1547(l0^°)s'^ 

+ .3035(10^^)3® + .4115(10^^)3® 

+ . 3672(10^® )s'^ + .1994(10^'^)s® 

Oil '2 9U 

+ .5782(10 )s t .7456(10 )s 

2U 

+ .2988(10 )/A(3) 

Gj^ 3 (s) = C- .1058(10®)s^® - .1135(l0®)s^'^ 

- .3974(10®)s^® - ,6970(l0^®)s^^ 
-•.7280(10^^)3^^ - .4933(10^'^)3^° 

- .2269(10^®)s^ - .7254(10^^)3® 

IQ 7 90 fi 

- .1622(10 )s - .2521(10 )s 

21 5 22-4 

- .2663(10 )s - .1835(10 )s 

- .7703(10^^)8® - .177l(l0^®)s^ 

- .1919(10^®)8 - .6876(10^^)]/A(s) 
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IS 114- 4 13 

G 2 ^(s) = C + 42.91s - .4241(10 )s 

. 2208(10"^ - . 8216(10® 

•fc .2254(10^^)3^® + .2870(10^®)s® 

± .1579(10^^)3®+ .5023(10^®)s’ 

■+ .1003(10^®)s® + .1285(10^®)s® 

on ii 91 

+ .1044(10 )s + .5152(10 )s 

+ .1408(10^^)3^ + .1784(10^^)3 

■ + .7423(10^®)']/A(s) 

622 ( 3 ) = C- 546.13^^ - .4005(l0®)s^^ - . 7l66(l0®)s^® 

- .523(10^®)s^^ - .1456(10^^)3^^ 

+■ .4209(10^®)s^® + .7393(10^®)s® 

+ . 5099(10^ ’^)s® + .21934(10^®)s^ 

+ .5989(10^®)s® + .1032(10^^)3® 

+ .1103(10^^)3*^ + .7002(l0^®)s® 

24 2 24 

+ .2399(10 - )s + .3529(10 )s 

24 

+ .1594(10 )3/A(s) 

623 ( 3 ) = C- .065753^®’- 4420.33^*^ - . 8*978(10®)s^® 

- .9572(l0®)s^^ - . 6854(10^® )s^^ 

12 -10 14 9 

- .3658(10 )s - .13769(10 )s 

- .3029(l0^®)s® - .1933(10^®)s'^ 

+ .-6354(10^''^ )s® + ..1416(10^®)s® 
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+ .7046(10^®)s^ - .5499(10^°)s^ 

- .5667(10^^)3^ - .1230(10^^)3 

- .6916(10^^)]/A(s) 

*l C T h IQ 

Gg^(s) = [ .0099s + 10.77s + 3779.8s 

+ . 6525 ( 10 ^) 3 ^^ + . 6614 ( 10 ^) 3 ^^ 

+ .4313(10^°)s^° + .1901(10^^)3® 

+ .5820(10^®)s® + ,12475(10^®)s^ 

+ .1860(10^^)3® + .1885(10^'^)s® 

+ .1244(l0^®)s^ + .4989(10^®)s® 

IQ 2 IQ , 

+ .1097(10 )s + .1134(10 )s 

+ .3845(10^®)]/A(s) 

632 ( 3 ) = [- 98.39s^® - .1032(10®)s^^ - . 3561(10® )s^® 

- .6111 (10^®) 3^^ - .6197(10^^)5^^ 

- .4051(10^^)3^° - .1789(10^®)s® 

- .5461(10^*^)3® - . 1161(10^® )s*^ 

- .1706(10®°)s® - .1691(10^^)3® 

00 ii 99 ? 

- .10817(10 )s - .4162(10 )s 

29 9 99 

- .8719(10 )s - .8677(10 )s 

- .2935(10^^)]/A(s) 

633 ( 3 ) = [ .50693^® + 532. 9s^^ + .1898(10®)s^® 

+ .3319(10®)s^^ t .3389(10^°)s^^ 

+ .2208(10^^)3^° + .9625(l0^°)s® 
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+ .2873,(l0^^)s® t .5909(l0^®)s'^ 

.+ .8299(iq^'^)s^ + .7754(10^®)s^ 

394 20 3 

+ .4591^(10 )s^ t .1607(l0^'')s^ 

20 2 20 

+ .3071(10 )s + .2841(10 )s 

19 

+ .9280(10 )]/ (s) 

where 

A(sO = + 1063. 8s,^^ + .3780(10^)3^^ 

+ .6691(10^)3^^ + .7021(l0^°)s^^ 

+ . 4777 ( 10 ^^) 3 ^^ + . 2215 ( 10 ^*^) 3 ^° 
+ .7195(0^^)3° 1 .1658(10^^)3° 

+ .27I5(l0^°)s^ + ,3125(10^°)s° 

20 5 21 4 

+ .2474(10 )s + .1297(10 )s 

+ .4257(10^^)3°' + .7998(10^^)3^ 

+- .7430(10^^)3 + .2411(10^*^) 


1. Two Input - Two Output Model 
For the Sixteenth order state model with the two inputs and two 
outputs indicated above, the transfer matrix becomes - 


G(s) -- 


Gii(s) • G^2 ^s) 


^2^(3) 622(5) 


Using the dominance algorithm in an INA mode over the frequency range 
0 ■_< (0 j< 200 with 200 equally spaced increments, system dominance was achieved for 
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.555 


0 


L = I 

0 1.0 

0 3840 

1.0 ■ ■ 0 

This result was obtained when the dominance control parameters -were set to 

unity. For 0^ = 0.0, dominance levels of [.74024, .96677] were obtained 

for rows 1 and 2, respectively, with 

^ .816 2289.0 

= 

1.0995 -47.39 

Since the Gershgorin circles provide no information apart from a. bound 

A 

on the eigenvalues of Q(s) (or Q(s) in a DNA), it is only necessary to 
compute the envelope of the Gershgorin and/or Ostrowski bands. For this 
purpose, the numerical method developed by Crossley [36] is used to calcu- 
late the envelopes centered about the corresponding diagonal element. 
Figures .4. 3 and. 4.4 display the Gershgorin bands for 
Q(^) = K^GCs)L 

over the frequency range 0 ^ w _< 200. 

Using precompensator the Gershgorin bands in Figures 4.5 and 4.6 
are obtained. These bands can be immediately reduced using the Ostrowski 
shrinking factors with zero feedback gains and are presented in Figures 
4.7 and 4.8. Note the high degree of band reduction between Figures 4.6 
and 4.8 which can be further reduced by increasing the feedback gain in the 
second loop. 

A 

To further examine the Ostrowski band for 922* frequency range was 
reduced to 0 £ oj ^ 50 arid is presented in Figure 4.9. From this figure it 
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Figure 4.6 Gershgorin Band for Row 2 using K 2 ^ .1 200 
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Figure 4.7 Ostrowski Band for Row 1 using K2 ^ 200 
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Figure 4.9 Ostrowski Band for Row 2 using K2 0 _< io _< 50 


is evident that the corresponding feedback gain is restricted to 
0 f 2 ^ 40- for stability in the. sense of Rosenbrock. 

Although ho further design attempts were made for the' two input case, 
it is apparent from the Ostrowski bands that dynamic compensation in the 
feedback loop may be required to achieve a step response with ■ minimal 
overshoot . 


2. Three Input t Three Output Model 
Using the G(s) matrix as 


GCs) = 


® 11 ^^^ 

G^iCs) G22 (s) G 23 (s> 

Sl^®^ 


an INA design with constant conpensators was initiated. When the post- 
compensator matrix was set to the identity matrix, diagonal dominance 
could not be obtained in all rows simultaneously. with the algorithm of 
section 3 . The next effort was to constrain, the 'L matrix to the form 


L 




0 0 


0 

0 



symbolizing the desire to control the individual outputs rather than linear 
combinations of the system outputs. With this structural form, row domi- 
nance was obtained in two passes through the basic- algorithm with the 
results 


0 .05 0 


L 


.01 0 0 


0 0 1.0 


and 
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1.000 1636.43 -1.3583 

K = .0001379 1.000 -.087418 

.02308 -206.685 1.000 

For this case, the dominance control parameters were selected as 
6^ = 0.75 0^ = 0.40 0g 0.10 

The co3?responding Gershgorin bands are contained in Figures 4.8 to 4.10 

I 

over the frequency range of 0 ^ m _< 50. 

Selecting the^ feedback parameters as 

150' 0 0 

F. ' = • • 0 5.0 0 

0 0 5.0 

the Ostrowski bands of Figures 4.10 to 4.12 were obtained. 

Using the K, L and F matrices indicated above, the interaction indices 
become 

= -0.2703 

I 2 = +0.0402 

I_ = -.00180 

o 

The first index suggests a moderate degree of system interaction in loop 1 
resulting from the closure of the remaining loops. This level of inter- 
action might have been anticipated from the level of dominance requested 
(0^ = 0.75). The negative sign for implies a constructive form of 
interaction in that closure of loops 2 and 3 augment the design efforts of 
loop 1 in the control of output y^(t). This design information can not be 
obtained from the -.Gershgorin or Ostrowski plots. 

The interaction index for output suggests a low level of inter- 

action upon closure of loops 1 and 3, The plus sign for l 2 » although not 
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Figure 4.10 Gershgorin Band for Row 1 0 50 
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Figure 4.15 Ostrowski Band for Row 3 0 , < < 50 



significant in this case, suggests that loop closure could have a deleter- 
ious effect on the control of y 2 (t). Thus, if the ma^itude of had 
suggested a higher level of system' interaction, design efforts to de- 
couple ^ 2 ^^^ could be considered. In this example, an alternate approach 
might be to restructure the post compensator matrix and re-examine dominance 
using the proposed algorithm. 

The interaction index for output y 3 (t) indicates a decoupled condi- 
tion. This result is in correspondence with the level of the dominance 
control parameter (9g = 0.1) specified and achieved by the dominance algo- 
rithm. 

To obtain the step response for the preliminary design above to com- 
manded step changes in the system outputs, the block diagram of Figure 
4.16 was used. 


/comm 



Figure 4.16 Closed loop system 
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In state variable form. Figure 4.16 becomes 
X = (A-BKFLC)X + BKF Y 
y{t) = C X(t) 

with Figures 4.17 to 4.19 representing the step responses corresponding to 
the following commanded input vector 

10. d 
10.0 

In each figure, the step responses have been superimposed to demonstrate 
that system interaction has been significantly reduced. 

To calculate steady state offsets in the step responses from the 

Ostrowski INA diagrams, Figure 2.4 of Section 2 may be used. From 

1 

Figures 4.13-4.15 

~150 0 o' 

F = 0 5 0 

_ 0 0 5 

the offsets are ' 


Bow 

1: 

100 = 


28.2% 

Row 

2: 


°;^(ioo)% = 

5.66% 

Row 

3; 



= 74% 


Thus the steady state values will be 71.8%, 94.34% and 26% of the commanded 
values. These values are easily verified from the step responses of 
Figures 4.17 to 4.19. 
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To demonstrate that the above design for the FlOO with three inputs 
and three outputs is not structurally unique, a new post.corapensator was 
selected as 

0 . .05 0 

•L = .0 0 2 

.01 0 0 

The initial guess for the precompensator was selected as the identity ma- 
trix with 0 .< m < 200 as the frequency range of interest subdivided into two 
hundred equally spaced points. 

Application of the algorithm provided diagonal dominance in all rows 

with 

.50275 

0 = " .22111 

.30606 

and 

2.9757 4834.2 -8.7416 

K = .001416 -13.906 .0696 

.005564 21.099 -1.5824 

The Gershgorin bands for 0 _5. “ 200 are contained in Figures 420 to 4.22 • 

Figures 4.23 to 4.25 show the same bands over the frequency range 0 < w < 50. 

Once diagonal dominance has been obtained, Ostrowski's theorem will 
apply for any set of 'stable feedback gains. Using the information that the 
FlOO is open loop stable, .all feedback gains were set to zero. The corre- 
sponding Ostrowski. bands are provided in Figures 4.26 to 4.28. From these 
figTXPes it is evident that the Gershgorin bands have been significantly 
reduced. Using the dominance levels above, the following minimum levels of 
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of reduction obtained are 


69% reduction of Gershgorin Band # 1 
50% reduction of Gershgorin Band # 2 
50% reduction of Gershgorin Band # 3 
The designer is now in a position wherein the Ostrowski bands can be 
further reduced by simply increasing the gains in each loop and recal- 
culating the shrinking factors for the selected set of gains. To complete 
the design, eigenvalue checks and/or step responses may be used. 

C. Discussion 

The applications considered in this section demonstrate the versa- 
tility and effectiveness of the dominance algorithm described in section 3. 
The algorithm is computationally fast and efficient with most applications 
requiring 100 CPU seconds or less to achieve the dominance condition. 

When the algorithm was tested against previously known results, two 
specific conditions were examined. First the reported results for the 
dominance producing compensators were implemented and verified for each 
case. In every instance, exact duplication of the gain space and Gershgorin 
and/or Ostrowski band was achieved. This condition established the accuracy 
of the algorithm in a non-optimization mode. The second condition ignored 
the reported parameter values and attempted to achieve diagonal dominance 
using alternate starting values for the compensators and the generalized 
optimization unit. Many new and interesting solutions were obtained and are 
reported in subsection A above. 

As a final test for the dominance algorithm, the sixteenth order state 
model for the FIDO turbofan engine was used to generate the appropriate 
transfer matrices corresponding to the two input and three input frequency 


101 



domain models. The results obtained clearly demonstrate the utility of 
the MNA design philosophy as a viable alternative for the' design of feed- 
baclc control units for the turbofan engine. Although the -results pre- . 
sented are preliminary, it is apparent that acceptable system performance 
using the design philosophy of section 3 is easily achieved. 
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SECTION 5 


CONCLUDING REMARKS 


The Multivariable Nyquist Array as defined in this report is the 
union of two mutually exclusive design methods; the Inverse Nyquist 
Array and the Direct Nyquist Array. The two design methods are mutually 
exclusive in -the sense that a design initiated in the inverse polar 
plane cannot in general be completed in the direct polar plane. This- 
apparent inconsistancy is due in part to the lack of duality between 
the definition of dominance for the INA and the corresponding definition 
for the DNA since 


''-1 

q. . f q. . 

However* the methods are structurally similar and thus provide the basis 
for the proposed design merger. 

Exploiting the' structural similarities between the two design methods, 
the dominance seeking algorithm of section 3 is appropriate for use in 
either the INA or the DNA design mode. In addition, constant or frequency 
dependent compensators can be evaluated for either row or column dominance 
with an indicated degree of preferred dominance. System interaction is 
easily assessed .and a means of compensator comparison has been provided. 

In its present form the MNA design algorithm performs in a batch com- 
puter mode. It is computationally efficient as demonstrated in section 4 

and provides an effective alternative design for turbofan engine control 
/ 

systems. In addition, the dominance algorithm' is ideally suited for imple- 
mentation on an interactive computer network. In this computer mode, it 
is conceivable that a complete design via the MNA could be accomplished 
within one working day. 
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